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£■ — . Abstract. We review the significant role played by red supergiants 

(RSGs) in stellar populations, and some challenges and questions they 

' j ' raise for theoretical stellar evolution. We present how metallicity and 

rotation modify the way stars go to the red part of the Hertzsprung- 
Russell diagram or come back from it, and how RSGs might keep a 
trace of their main-sequence evolution. We compare theoretical popu- 
O | lation ratios with observed ones. 
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1 Introduction 



1.1 Why red supergiants are interesting objects 



If we make the hypothesis that all stars with a mass between 8 and 25 Mq go 
through a red-supergiant (RSG) phase, it means that 80% of all massive stars 
will be a RSG one day. This phase is a crucial one for the mass-loss history of 
the star (see C. Georgy 's contribution in this conference and references therein). 
t-h RSG are also supposed to be the direct progenitors of the most numerous type of 
core-collapse supernovae (cc-SNe): the type II. 

We present an evolutionary scheme for massive stars in Table [TJ inspired from 



C*~) Conti (1975), based on solar-metallicity non-rotating models. In this scheme, all 
1— I massive stars with M < 30Af Q end their life as a RSG, while all massive stars 
*^ with M > 40 M Q end their life as Wolf-Rayet (WR) stars. Between 30 and 40 M , 
r |~j the mass loss during the RSG phase might be sufficient to strip the external layers 
rS deeply enough for the star to become a WR star. It must be clearly stated here 
that the real evolutionary path between the different stellar types depends on 
the physics considered in the models, and might be significantly modified by the 
mass-loss rates used, the overshooting amplitude, the inclusion or not of the effects 
of rotation and of magnetic fields. A large fraction of massive stars might also 
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M > 60 M : 

M = 4O-6OM : 


-> Of/WNL -j- LBV -> WNL -> (WNE) -> WC 
-> BSG -> LBV -> WNL -> (WNE) -> WC 


-}• SN Ibc 
-> SN Ibc 


WR 


M = 


3O-4OM : 


-> BSG -> RSG -> WNE-s- WCE 


-> SN Ibc 




M = 
M = 


25-3OM : 
1O-25M : 


-> (BSG) -)- RSG -> (YSG?) 

RSG -> (Ceph. loop for M < 15 M ) RSG 


-> SN II-L/b 
-> SN II-P 


RSG 



Table 1. Modified Conti scenario for the evolutionary scheme of massive stars. 



undergo interactions in close binaries (Sana et al., 2012), and thus follow a quite 
different evolutionary path. 



RSGs are easily seen in external galaxies (Langon et al., 2009 Drout et 



2012 Neugent et al., 2012) which makes them also interesting objects because 



they can be used as metallicity or distance indicators. Davies et al. (2010) and 



Bergemann et al. (2012) proposed the RSGs as metallicity indicators, showing that 
the determination of Z is robust even when only modest resolution (R — 1 — 3000) 
is available. They are proposed as distance indicators (see B. Davies' contribution 
in this proceedings). Since the RSG phase is limited in mass range and is a 
relatively short advanced stage of stellar evolution (< 10% of the lifetime), they 



are also proposed as age indicators (Langon et al., 2009) 



1.2 Some present-day challenges in the physics of RSGs 

In ID stellar-evolution codes, many mechanisms cannot be derived from ah initio 
physical principles, and are thus parametrised. It is the case for a mechanism that 
is of extreme importance for RSG modelling: convection. 

Looking at the structure of a RSG, we see that while the core is encompassed 
in roughly half a solar radius, the star itself might extend to around 800 Rq, most 
of it being convective. The modelling of such objects makes it necessary to have 
a good treatment of convection, which is still out of reach of ID stellar-evolution 
codes. This convection is non adiabatic, and supersonic turbulence needs to be 



taken into account (Maeder, 2009) 



Another challenge concerns the mass loss. It is implemented in the codes thanks 
to prescriptions offered in the literature. However there is a debate about the rates 



inferred for RSGs. Some observations (Mauron & Josselin, 20111 give support to 



the de Jager et al. ( 1988 1 rate usually implemented in the codes for this stage of 
the evolution. Others, on the contrary, challenge these rates and suggest much 



larger ones (van Loon et al., 2005). A large uncertainty comes from the unknown 



nature of the mass loss in that stage: steady and regular, or outbursting? If the 



mass loss occurs by bursts (Humphreys et al, 2005 Smith et al., 2009 Moriya et 



al., 2011), the probability of observing a star precisely in a burst episode might 
be quite low, so the mass-loss rates are probably underestimated. In the stellar 
modelling perspective, we are bound to use steady rates and thus implement a 



RSG evolution 



3 



mass-loss recipe that would average over time the total mass-loss episodes gone 
through by the RSG. However it is unclear how to derive this time-averaged mass 
loss out of burst episodes, and moreover the steady approximation might modify 
the global behaviour of the star. 



1.3 Some crucial questions for the evolutionists 

The evolution to or from a RSG stage is not straightforward to understand. We 
need to determine the physical conditions driving a star to be or not to be a RSG 
through the following questions: 

• When does a massive star enter into the RSG phase? 



• How long does this phase last? 



• What is the impact of the physical processes available in stellar evolution 
codes (treatment of convection, mass loss, rotation, magnetic fields, multi- 
plicity)? 



• What types of cc-SNe arise from RSG? 



- type II-P? II-L? 

— Is there a SN event if the core collapse results in a black hole? 

In Section [2j we shall address the first three questions by studying theoretical 
models published recently (Ekstrom et al, 2012 Georgy et al., 2013). In Section[3j 
we shall consider RSGs as a population and compare theoretical ratios to the 
observed ones. 



2 To be or not to be a RSG 



Let us first quickly recall the main physical ingredients of the models (the details 
can be found in Ekstrom et al. 2012). Compared to the previous complete grids 
of Schaller et al. (1992), some reaction rates have been updated, the abundances 
have been chosen to be those of Asplund et al. ( 2005 1 , and the opacities have 
been updated to match the abundances. The effect of rotation are implemented 
according to Zahn ( 1992 1 for the horizontal turbulence diffusion coefficient and to 
Maeder (19971 for the shear diffusion coefficient. 

Mass loss is an important ingredient for the modelling of massive stars. In the 
models, we use Vink et al. ( |2001[ ) for the O-type stars and for the RSGs we use 
a linear fit based on observations by Crowther (2001) that yields mass-loss rates 
similar to the de Jager et al. (1988) prescription. For stars with M > 20Af Q , we 
multiply the mass-loss rate by a factor of 3 if the star has some external layers 
that have a supra-Eddington luminosity. 
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central 4 He mass. frac. central 4 He mass. frac. 

Fig. 1. Evolution of the T e g as a function of the central He combustion (time evolution 
goes from left to right) in non-rotating models. The RSG domain is shaded in red. The 
models of 9Mq (purple), 12 (dark blue), 15 (cyan), 20 (green), 25 (yellow), 32 (orange), 
and 40 Mq (red) are represented. Left: solar metallicity (Zq = 0.014). Right: SMC 
metallicity (Z = 0.002), models from Georgy et al. (in prep). 

2. 1 Becoming a RSG 

At solar metallicity (Z = 0.014), our models cross quite rapidly the Hertzsprung- 
Russell (HR) gap, spending most of the core He-burning phase in the red-supergiant 
region (see G. Meynet's contribution in this conference for a discussion on the 
physical mechanisms involved). 

Figure [T] shows the behaviour of the T e g as a function of the central He 
abundance for two different values of Z. The mass range for becoming a RSG 
(log(T c ff) < 3.70, red-shaded area on the figure) is between 9 and 40 M for our 
non-rotating models at Z (left panel), and the associated age range is between 
4.5 and 30 Myr. It appears that the 9 and 12 M models cross the HR diagram 
right after the main sequence (MS), before He ignition. The 15 and 20 M models 
cross the gap before having burnt half of the central He and finish their life in the 
red. The 25 M model crosses the gap slightly earlier but then makes two blue 
loops during core He burning, and finishes its life in the blue. Both the 32 and 
40 Mq models make a quick incursion into the red, then make a large blue loop, 
come back very briefly and finish their life as WNE. 

As discussed in G. Meynet's contribution in this conference, when the metal- 
licity is lower, the star enters into the RSG stage at a more advanced stage of the 
core Helium-burning phase (Fig. [I] right). This implies that the duration of the 
RSG phase is much shorter. While at solar metallicity, RSGs are He-burning stars, 
at low Z they are C-burning stars (see J. Groh's contribution in this conference). 
We shall see below the consequences it bears on stellar populations. 

Rotation tends to keep the most massive stars in the blue, so the mass range 
in which we expect to see stars evolve into RSGs is reduced in our V/Krit = 0.40 
(fi/f2 cr it = 0.56) models: between 9 and 25 M only. The age range associated 
with this phase is between 8 and 35 Myr. Looking more closely at the effects of 
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Fig. 2. Same as Fig. [Tj for 15 Mq models. Left: solar metallicity [Z = 0.014), various 
rotation rates (models from Georgy et al. 2013). Right: SMC metallicity (Z = 0.002), 
various rotation prescriptions (models from Meynet et al. 2013, with fi/fi C rit = 0.50 
initially) . 



rotation, we note that the trend changes depending on the rotation rate. Passing 
from VyV^rit = to 0.20 (fi/fi C rit = to 0.3) favours the redward crossing of the 
HR diagram. But if rotation rates become higher, the trend gets inversed and the 
crossing occurs at later times during He burning (Fig.|2j left, models from Georgy et 
al. 2013). In real stellar populations, we expect a velocity distribution, and thus a 
diversity of behaviours. One must bear in mind that there are several prescriptions 
for implementing the effects of rotation into the codes. In the literature, we find two 
different shear-diffusion coefficients (Maeder, 1997 Talon & Zahn, 1997) and three 



different horizontal-turbulence coefficients (Zahn, 1992 Maeder, 2003; Mathis et 



\al, 2004] ) . Depending on the set of prescriptions used, the behaviour might be 
quite different (Meynet et al, 2013), as illustrated in Fig. (right). 



2.2 Back to the blue 



The duration of the RSG phase depends strongly on the mass loss experienced by 
the star (see C. Georgy's and G. Meynet's contributions, and references therein). 
According to the observations of Smartt et al. ( 2009 ) , there are no RSG progenitors 
for SN II-P above around 18 M© . This coincides with the increased sensitivity to 
supersonic convection at 20 M and above (Maeder, 20091, and thus potentially 
higher mass- loss rates. In our rotating models at Z@, the lower-mass stars end 
their lifes as RSGs and the higher-mass stars move back to the blue. The transition 
between both behaviours occurs near 18 Mq. In the models, it is most probably 
due to the treatment of the supra-Eddington layers, however it illustrates well the 
role played by the mass loss. 
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Fig. 3. 15 M© models at Z@, from Georgy et al. 2013. Left: evolution of the surface 
velocity as a function of the T e s during the HR-gap crossing. Right: N/C surface abun- 
dances (in number) as a function of the surface velocity in the RSG phase (same colour 
coding as Fig. [2] left). 



2.3 Keeping a trace from the past 

One can wonder whether a rotating star might keep a trace of its former rotation 
once it has undergone the tremendous inflation leading to the RSG phase. If the 
crossing of the HR gap is accomplished within a Kelvin-Hclmholtz timescale (which 
is the case for models crossing right after the MS and before any He burning), the 
surface velocity of the star will evolve according to angular-momentum conserva- 
tion. On the contrary, if the crossing occurs during part of the He-burning phase 
(which is the case for the most rapid rotators), the winds have time to remove 
some angular momentum during the crossing, but also the transport mechanisms 
have time to bring some angular momentum to the surface, and the evolution of 
the surface velocity is more complex (Fig. [3j left). 

Once in the red, we could expect the subsequent dredge-up to erase any previ- 
ous enrichment, but this is not the case. Figure [3] (right) shows that the enrichment 
of the surface of a former rapidly-rotating star remains more highly enriched (by 
a factor of 50-60) than that of a non-rotating star. We can distinguish three end- 
point zones: one were stars have very low surface velocity (V cq < 0.1 km/s) and 
surface enrichment (N/C < 7, green zone in the figure); one with surface velocities 
around 0.2 km/s and N/C around 10 (blue zone); and the zone where the former 
most rapidly-rotating stars end, with velocities around 0.5 km/s and N/C > 15 
(red zone). 



3 RSGs as a population 

Stellar-population ratios are good indicators to compare theoretical tracks with 
"real" stars. There is a difference between ratios observed in a coeval population 
(open cluster) or in a general population with a constant star-formation rate: 
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log(age) time [Myr] 



Fig. 4. Left: Evolution of the population ratio between red evolved stars and MS stars (in 
an interval of 2 mag below the turn-off) (theoretical curve: black solid line; observations 
in open clusters: grey squares, size and filling according to the MS star sampling). Right: 
Theoretical evolution of the number of WR (black line) and RSG (grey line) in a cluster 
of 10,000 stars at birth. 

• in a coeval population, the ratio between two types of stellar populations 
(A and B) can be predicted by: 

B ImbZ™ ^( m ) dm 

where we just have to know the initial mass range leading to each population 
[M Aymin ...M A . max \ and [M B , m i n —M B , ma ^\ and integrate it over an IMF 
(j>(m). It thus brings constraints on the mass limits inferred; 

• in a constant star-formation population, the same ratio involves also 
the duration of the phase in which the star is of the given type (ta and tb): 

It thus brings constraints on the mass limits and the durations. 
3.1 RSGs in coeval clusters 

The evolution of a population of red evolved stars in a cluster can be understood 
in the view of stellar evolution as follows (Fig. |4j left): 

• at very young ages, no star has yet had time to become a RSG, so none is 
expected; 
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• around 10 Myr, one expects the number of RSGs to rise steeply, with a peak 
when high-L RSGs (i.e. those arising from the most massive models going 
through a RSG phase) dominate the population; 

• a drop-down occurs in the age range of wide-loop Cepheids; 

• the number increases again when low-mass red giants start dominating. 

In Fig. [4] (left), we evolved an initial population of 10,000 stars and counted 
as RSG the Mj n i > 9M stars having log(T e ff) < 3.70 and as MS stars all non- 
evolved stars populating an interval of 2 magnitudes below the turn-off (TO). 
When we compare the theoretical curve to observed open clusters, we expect 
a large dispersion because we run into a stochasticity induced by small-number 
statistics. Typically, all clusters with a point falling well below the theoretical 
expectation contain less than 20 MS stars (2 magnitudes below the TO), the only 
exception being the NGC 3532 + NGC 6494 point. 

Figure [4] (right) shows the evolution of WR-stars and RSGs populations respec- 
tively. It is interesting to note that on theoretical grounds, we expect practically 
no overlap between both populations in the frame of single-stars evolution. Hence 
if some WRs would be observed in a cluster of RSG age (between around 10 to 20 
Myr), it could provide a test for the binarity channel associated to the formation 



scenario of WR stars (Eldridge et al., 2008[ and references therein). The opposite 



case (RSGs in clusters of WR age, i.e., younger than 8 Myr) would rather chal- 
lenge the homogeneity of the population (cluster is not coeval, chance alignment, 
...). 

3.2 RSGs in constant star-formation regions 

The study of general populations is often done throughout a range of metallicities, 
and thus of different galaxies. In this case, the RSGs counts are performed in a 
restricted luminosity range (log(L/L Q ) > 4.9), in order to avoid a contamination 



by red giants. We thus adopted the same criterion, using the Ekstrom et al. ( 2012 ) 
grids and similar grids at Z = 0.002 for which publication is in preparation. 

Massey ( |2002[ ) reports the values of the RSG/WR ratio observed in the metal- 
licity range going from M31 (12 + log(0/H) = 9.0) down to the SMC (12 + 
log(0/H) = 8.13). He finds a linear increase of this ratio with lowering the metal- 
licity (see Fig. [5j left). When we compare these results to models, there is a 
problem. Without rotation, the predicted ratios are much too high (by 0.5 to 1 
dex), but the metallicity trend is well reproduced. With rotation on the contrary, 
the solar value is in relative agreement with the observations, but the metallicity 
trend is not reproduced and the SMC value is much too low. Note that we used 
the same initial velocity (V/V CT n = 0.40) as average rotation rate than for the Zq 
grids, and it might be too rapid for the average population at Z — 0.002. A slower 
rotation rate might better reproduce both the values and the Z-trend. 



Looking at the BSG/RSG ratio as reported by Eggenberger et al. ( 2002[ ) (Fig.[5| 



right), the situation is not better. While at solar metallicity, the values of both 
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Fig. 5. Population-ratio predictions from the Ekstrom et al. (2012) models (grey pen- 
tagons: non-rotating, open symbol; rotating, filled symbol) overplotted on observational 
results from the literature. Note that the figures are exactly as published: when a theo- 
retical point falls beyond the boundary of the plot, it is replaced by an arrow. Left: RSG 
to WR ratio. Figure 11 of Massey (2002). Right: blue- to red-supergiant ratio. Figure 1 
of Eggenberger et al. (2002, black triangles: 'B' contains 0-, B-, and A-type stars; black 
circles: 'B' contains only B-type stars). 



rotating and non-rotating models are quite in agreement with the observations, 
we definitively lack RSGs at low metallicity. The values for the SMC are much 
too high (94 for the non-rotating models, and 87 for the rotating ones, beyond the 
figure's boundary). This reflects the late crossing of the HR diagram at low Z. 
Note that we saw in Sect. |2.l| that varying the rotation rate had various outcomes 
on the time of the crossing. Based on duration considerations, we find for the 
15 Mq rotating at V/V crit = 0.40 a BSG/RSG rati o of 6 5, while if we consider 
a velocity distribution as proposed by Huang et al. (2010), the same ratio drops 
down to 16, still too high but not as dramatically different. It points nonetheless 
to a possible solution to the population-ratio problem: by having the stars crossing 
the HR gap much earlier after the MS, which would increase the duration of the 
RSG phase. 



4 Conclusion 



The RSG phase is a delicate one for theoretical stellar models, because stars in 
this phase are dominated by mechanisms that are a challenge for ID computa- 
tions: convection and mass loss. Improvements in the observations of the RSGs 
mass-loss rates are highly needed to get reliable prescriptions to include in the evo- 
lution codes. Developments brought by 3D modelling of convection are becoming 
available and might bring physical prescriptions that can be included in the ID 
codes. 
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The conditions for a star to become a RSG are sensitive to the structure at the 
end of the MS, mainly the size of the core and the amplitude of its contraction 
after H exhaustion. It is modified by overshoot, rotation, metallicity. Comparison 
of stellar-population ratios brings an observational constraint that points to a lack 
of RSGs produced at low metallicity by the models. 

Let us conclude by quoting Langer & Maeder ( |1995[ ), who finish their paper en- 
titled precisely The problem of the blue-to-red supergiant ratio in galaxies 
with these words: "We have to conclude that massive star models in the consid- 
ered mass range still lack some significant physical ingredient. However we want 
to emphasise that this does not imply that the results of massive star theory have 
to be questioned altogether. The B/R-ratio is a quantity which is known to depend 
extremely sensitive on the model parameters. It is thus a welcome amplifier which 
can (and finally will) be very useful to constrain the model physics very accurately." 

We see that nearly 20 years later, the improvements brought to the models 
are not yet able to solve this longstanding problem. On the point of view of the 
physics, we understand the behaviour of our models, but Nature seems not to be 
eager to follow these rules. Since Nature is always right, we need to explore further 
the challenge she addresses. 
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